Spin-polaron excitations in a doped Shastry-Sutherland model 
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Using variational algorithm on an infinite Shastry-Sutherland (SS) lattice we show that the in- 
troduction of a static nonmagnetic impurity into a dimerized ground state leads to a formation of a 
small, localized spin polaron surrounding the impurity site. Due to a particular symmetry of the SS 
lattice, the polaron is extremely anisotropic with a short spatial extent. The presence of nonmagnetic 
impurities leads to a formation of pronounced in-gap peaks in the dynamical spin structure factor, 
which we attribute to the spin-doublet excitations of a single unpaired spin S — 1/2 surrounded by 
triplet fluctuations. Our results are relevant for the description of SrCu2(BOa)2 compound when 
doped with nonmagnetic atoms at Cu sites. 

PACS numbers: 75.10.Jm, 75.30.Hx, 75.40.Gb, 75.50.Mm 
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I. INTRODUCTION 

The relation between disordered spin liquids with a 
gap in the spin excitation spectrum and superconductiv- 
ity has aroused a lot of interest. It was suggested^ that 
doping a gapped system may lead to hole pairing and su- 
perconductivity. In this sense the idea of finding a super- 
conducting state in a doped two-dimensional spin-liquid 
SrCu2 (1303)2 compound, representing the first physical 
realization of the Shastry-Sutherland (SS) modelf 2 ^^ has 
been proposed^*! due to its structural similarity with 
the high-temperature cuprates. So far, numerical calcu- 
lations do not seem to support superconductivity in the 
doped SS model £ 

Another aspect of the study of the doped SrCu2 (1303)2 
is the influence of impurities on the nature of the singlet 
(dimerized) ground state. When J'/ J of the underlying 
SS model is close to but below the critical value ~ 0.7, 4 
at finite doping, the system may select other ordered or 
disordered states, which may in turn strongly affect the 
magnetism of this material. It is rather straightforward 
to realize that replacing magnetic Cu ions with nonmag- 
netic ones breaks the singlet nature of the bonds, thus 
providing an additional frustration, resulting in the re- 
duction of the spin gap. The effects of doping should 
materialize in inelastic neutron scattering through an en- 
hanced lifetimes of single-triplet excitations - magnons. 
A finite lifetime associated with magnon scattering on 
spin excitations surrounding the nonmagnetic ions would 
in turn lead to a broadening of the single-triplet peak line 
in the doped SrCu2 (603)2- 

The SS model consists of a relatively simple spin (S = 
1/2) Hamiltonian, 



H — J S Si • Sj + J' S Sj • Sj] 



(1) 



here (i,j) denotes the nearest neighbors (NN) on bonds 
J (dimers) and next-nearest neighbors (NNN) on 

bonds J'. We introduce a single nonmagnetic impurity - 
a site with a missing spin, into the model in Eq. l|T]l. It 



is well known that for J'/ J < 0.7 the exact ground state 
can be written as a product of singlet dimers: 4 We should 
stress that the introduction of a single hole into the exact 
dimerized ground state has nontrivial consequences. 

Previous theoretical studies have been focused mostly 
on the doping of SS model with mobile holes with 
the intention to explore the hole-pairing mechanism rel- 
evant for superconductivity. In our work we limit our 
investigations to one static nonmagnetic impurity (mag- 
netic polaron) within the SS lattice. Our calculations 
are relevant to the case where the substitution with non- 
magnetic ions takes place at Cu sites. We use a varia- 
tional method within a variational space constructed by 
successively applying the off-diagonal parts of Hamilto- 
nian on the starting (Oth-generation) state vector. 
The method yields converged results, valid within the in- 
finite system, for static quantities such as the value of 
the spin gap as well as for the spin distribution func- 
tion. We also explore the properties of the first few low- 
lying excited spin-polaron states. Finally, we compute 
the dynamical spin structure factor as a function of the 
wavevector. 

For simplicity, we have decided to neglect the 
anisotropic Dzyaloshinsky-Moriya interactions originat- 
ing in spin-orbit coupling 10 and used the (bare) SS 
Hamiltonian QJ. 



II. VARIATIONAL APPROACH IN 1-HOLE 
LIMIT 

A variational space was constructed by starting from 
the Oth-generation state vector |1q) representing one 
doped hole with a neighboring uncoupled spin within a 
sea of singlet dimers (Fig.^a)). First generation, M = 1, 
of variational states was consequently obtained by acting 
with H (JIJ on state 1 1 q ) . Successive generations leading 
to Mth order states {|«m}} were obtained starting from 
all previously generated ones {|«m-i)}. A state can be 
uniquely determined by the orientation of the spin adja- 
cent the impurity site, S* = ±1/2, and the arrangement 
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Figure 1: (color online) Two initial states (M = 0) in the 
1-hole variational procedure defined on infinite SS lattice: (a) 
state |lo) with uncoupled spin = 1/2 (|) within a sea 
of singlet dimers (black solid lines), and (b) state | 2q ) with 
S* = -1/2 (|) and one S z = 1 triplet dimer (TT)- Dashed 
black lines denote off-diagonal terms (J') of H 



of S z = 0, +1 and —1 triplet dimers in the vicinity of 
the hole, {r°}, {r^ 1 } and {r" 1 }, respectively. With rf* 
we denote the positions of triplet-dimer centers relative 
to the impurity site, i-th variational state within Mth 
generation is thus written as 



^I^.K},^ 1 },^,- 1 }) 



(2) 



where a fixed ordering of rf vectors is assumed to pre- 
vent double counting. In this notation the initial state is 
written as | lo) = 1 1/2, {}, {}, {}}. In a similar fashion we 
can also describe the (exact) 0-hole ground state. 

An advantage of the introduced procedure, in which 
a variational space grows like N st ~ 6 A/ (see Tab. QJ, is 
that triplet dimers of newly generated variational states 
develop gradually around a hole, which is very promising 
in a way that only a few such states would suffice to 
correctly describe the localized nature of the spin polaron 
found in Refs&S As application of the off-diagonal part 
(J') of Hamiltonian H on intcr-singlet-dimer bonds is 
zero, the relevant region from where new triplets may be 
generated is limited to a neighborhood of the hole. Due 
to specific topology of the SS lattice, a gradual expansion 
of triplet region with increasing M is proceeded in a zig- 
zag path. A short calculation shows that already M = 3 
is enough for the outermost triplet to exceed the square 
region of 32 sites used in Ref. 8 . 

Along with the construction of the variational space we 
calculated nonzero matrix elements Hfj = {iM\H\jM-i) 
separately for each term in Hamiltonian (JTJ . Matrix 
was diagonalized at each generation M to control the 
convergence. In Fig. [21 we show the convergence of a few 
lowest eigenenergies for the 1-hole and, for comparison, 
also for the 0-hole case. Concentrating first on the single 
impurity case, we plot the energy difference between the 
ground state and the first two excitations, and A^ 1 ), 
as a function of J' J J at fixed J = 76.8 K. In this par- 
ticular calculation two initial state vectors were used to 
improve the sampling of the variational space: | lo) and 
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Table I: Variational calculation for J = 76.8 K, J'/J = 0.62. 
M (N a t) denotes the number of variational steps (states) and 
A (0) = e[ 0) (S = 1) - E£ } (S = 0) represents the spin gap of 
the undoped system. Eq is the ground state energy of the 
1-hole system measured relative to the ground state energy 
E^ 0) (S = 0) of the 0-hole system. A (1) = e[ 1} (S = 1/2) - 
E^ iS = 1/2) and A« = E?\S = 1/2) - e£\S = 1/2) 
where e[ 1] (S = 1/2) and E^ ] (S = 1/2) represent the first 
and second excited states of the 1-hole system, respectively. 
We also show extrapolated values for M — > oo. All energies 
are listed in meV. 



|2 ) = |—l/2, {}, {(0, 1)}, {}> with one S z = 1 triplet 
located above the uncoupled spin as also depicted in 
Fig^b). In addition, we choose to work in the S z = 1/2 
sector that contains the ground state of the 1-hole system 
in the absence of the magnetic field. No other symme- 
tries were explicitly employed in the construction of the 
variational space. The SS model doped with a static hole 
possesses the mirror symmetry over the plane containing 
a doped dimer. We have checked that eigenstates and 
corresponding correlation functions transform according 
to this symmetry. 

Moving back to the convergence of results for the 1- 
hole presented in Fig. [3 we notice, that values of 

A' 1 ' and A^ 1 ) have nearly converged for J'/J < 0.62 for 
M = 8, i.e., their values at M = 8 are not far from 
their extrapolated values for M — > oo. As expected, 
the variational approach becomes worse when J'/J ap- 
proaches the transition to the AFM state of the 0-holc 
case, J'/J ~ 0.7. This can be better seen in a 0-hole 
case where the spin gap A(°) should close as J'/J -> 0.7. 
In our calculations J'/J ~ 0.64 represents a critical ra- 
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Figure 2: (color online) Energies relative to the ground-state 
energy vs J'/J calculated at J = 76.8 K for 0- and 1-hole 
system, defined on infinite SS lattice, using M steps of the 
variational method. Black solid line represents experimental 
value obtained on SrCu2(BC>3)2. In the inset we show the 
comparison with the 6th, 8th, 15th and 21st order dimer- 
expansion calculation in J'/J for holes 



tio beyond which the convergence of the spin gap is 
lost. On the other hand, at J'/J = 0.62 a single-triplet 
energy is quite well converged for M = 8, which in turn 
leads to thermodynamically converged value of the spin 
gap A(°) that also agrees with the experimental value 
A'°J ~ 3 meV. This suggests that the same set of J and 
J' as in calculations on 20 sitesi&iiiii may be used here 
in order to compare with realistic SrCu2 (603)2 system. 
We attribute this fast convergence to the thermodynamic 
limit to the fact that single-triplet excitations are well 
localized, which renders a limited variational space suf- 
ficient for nearly exact description of low-energy excita- 
tions. 

A 0-hole case may also serve as a test for the method, 
for it can be directly compared to alternative approxi- 
mate approaches. A comparison with dimer-expansion 
technique of Zheng et al}§ and Knetter et alH is shown 
in the inset of Fig. [3 Our method for M = 8 gives sim- 
ilar energy-gap dependence as the 21st order calculation 
for J'/J < 0.64, however, at larger values of J'/J > 0.64 
our method captures better the closing of the spin gap. 
It should be noted at this point that even better conver- 
gence with M is expected when a translational symmetry 
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Figure 3: (color online) Spin distribution function S„(r) cal- 
culated in (a) ground state, (b) first excited state, and (c) 
second excited state of the 1-hole SS model Q. Calculation 
was performed for J = 76.8 K, J'/J = 0.62 and M = 8. 



of the lattice is employed for the 0-hole case. 



III. RESULTS 

We have calculated the shape of a single spin polaron, 
forming in the vicinity of a static hole, as measured by 
the spin distribution function calculated in the 1-hole (by 
hole we denote the magnetic impurity) eigenstate In' 1 )) 
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Figure 4: (color online) g^-averaged 1-hole dynamical spin 

structure factor Siz (q x ,u) as a function of M calculated for 
J = 76.8 K and J' /J = 0.62. Lines are shifted for clarity. 
The arrows denote extrapolated positions of the lowest two 
peaks. 
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of H (HJ: 



S'{T) = (nV>\S' v \nV 



n = 0,l,2. 



(3) 



Here r denotes the position of the spin with respect to 
the hole. We first note that in the case of J'/ J < 0.7 and 
zero doping S§ (r) = for all r due to a singlet nature of 
the ground state |0 ( -°- ) ). 

In Fig. |3 we present the shape of the spin polaron, 
Sq(t), and of its two lowest excitations, £f(r) and 5*|(r), 
calculated for M = 8. The zero value of S„(r) indicates 
that r points either towards the edges of a singlet or a 
S z = triplet dimer. Note also that the following sum- 
rule holds: J2r^n( r ) = 1/2- From Fig. El it is apparent 
that the ground state of the spin polaron is rather well 
localized around the impurity. Only first two neighbors 
lying in the direction perpendicular to the dimer contain- 
ing the impurity are strongly affected by the presence of 
the impurity. The other two neighbors lying along the 
direction of the dimer with the impurity are nearly un- 
affected giving iSo(r) ~ 0- The shape of the spin polaron 
is therefore extremely anisotropic. This is in agreement 
with the exact T = calculations on 32 sites£ Accord- 
ing to Ref£, a mobile hole introduced to the dimerized 
ground state forces the resulting free spin to minimize 
the energy by forming a five-spin chain with two of its 
nearest dimers. This excitation is localized and has a 
small dispersion compared to the t-J model on a square 
lattice. These results are in qualitative agreement with 
our small-spin-polaron picture around the localized hole. 

In Fig. |3| we also show a spatial distribution of the 
lowest two excited states of the spin polaron, Sf(r) and 
S'l(r), calculated as well in S z = 1/2 sector and for 
M = 8. In contrast to the Sg(r) case, Sf(r) shows 
spin disturbance predominantly along the direction of 
doped dimer. Intriguingly, this excitation is not centered 



Figure 5: (color online) Intensity plot of the 1-hole dynami- 
cal spin structure factor Si* (q,L)) calculated for J = 76.8 K, 
J' I J = 0.62 and M = 8. A map was obtained by interpo- 
lation between 40 equidistant values of q x 6 [0,4]. Note also 
the logarithmic intensity scale. 



around the impurity; instead it is shifted away from it to- 
wards the neighboring orthogonal dimer. We note, how- 
ever, that this excited state of the polaron is not confined 
by the extent of the variational space which allows dis- 
placements at much larger distances. This leads to the 
conclusion that the lowest polaron excitation is centered 
on the neighboring orthogonal dimer, yet it remains lo- 
calized in the vicinity of the hole. A slightly different 
picture is obtained for S§(r) where the uncompensated 
spin is distributed more evenly on NN and NNN around 
the hole. As naively expected, higher energy excitations 
of the spin polaron spread further away from the impurity 
site. For this reason more variational states are needed 
to obtain full convergence in this case. 

In the last part we investigate the effect of spin polaron 
states, forming around static nonmagnetic impurities, on 
scattering experiments. In Fig. 01 we present the dynami- 
cal spin structure factor calculated within our variational 
approach at finite wavevector q, 



Sgfaw) = -Im/o* 1 



5= 



1 



S* 



o' 1 



Sq = E eiq(R+r) ^ 
R,r 



(4) 



Here R runs over all unit cells reached within M steps 
and r spans four vectors forming the basis of the unit 
cell that contains two orthogonal dimers. For the details 
describing interatomic distances in SrCu2 (603)2 system 



5 



we refer the reader to Refii&. In addition, a finite e = 
0.25 meV (cj + = oj + ie) is used to smooth the spectra. 

In Fig. 0] we show spectra of SzJ(q x ,w) averaged over 
the interval of q x S [0,4] (in units of 27r/ao) for vari- 
ous M. It is evident that spectra at w < 2 meV has 
fully converged, while convergence at higher frequencies 
is slower. Nevertheless, from Fig. 0] it is apparent that 
finite doping leads to a finite response of the dynamical 
spin structure factor at frequencies below the gap value, 
lj w 3 meV, of the undoped system. Similar effect is as 
well seen in the surface plot presented in Fig. [3] where 
we explicitly provide the q x dependence that was com- 
puted using 40 discrete values evenly spaced on the in- 
terval q x G [0,4]. Despite slower convergence at higher 
frequencies, as observed in Fig. 01 two predictions can 
be made based on presented calculations. We first pre- 
dict that in-gap peaks should appear at frequencies far 
below the gap value in the scattering experiment on a 
doped SrCu2 (1303)2 compound, where some Cu atoms 
would be replaced with static nonmagnetic impurities 
(holes). These in-gap peaks would carry information 
about the excited spin-polaron states. Second, doped 
impurities would contribute to broadening of the one- 
magnon peak at u> ~ 3 meV through inelastic scattering 
of magnons on excited states of spin polarons surround- 
ing the impurities. This effect can be clearly seen from 
Figs. 0] and [S] where extra peaks appear in the vicinity 
of u> ss 3 meV in comparison with calculations on the 
undoped system i 15 ' 17 



IV. CONCLUSION 

In summary, we presented static and dynamic proper- 
ties of a single spin polaron in a doped SS lattice. We 
have developed a variational approach, defined on the 
infinite lattice, that gives numerically highly accurate re- 
sults in the thermodynamic limit of static properties of 
spin polaron and its lowest excitations. The first and 
second excited state of spin polaron represent 5=1/2 
excitations, their energies lie well within the spin gap 
of the undoped system. The shape of the spin polaron 
is highly anisotropic while its extent is well localized in 
agreement with exact diagonalizations on finite systems. 8 
Spin polaron extends predominantly to the first neighbor- 
ing dimers lying in perpendicular direction with respect 
to the dimer containing the impurity. The shape of the 
first excited state of the polaron extends further in space 
and is off-centered with respect to the position of the 
impurity. Spin structure factor of the doped SS model 
is consistent with the appearance of in-gap peaks and 
the broadening of the one-magnon peak due to inelastic 
scattering of the magnon on static spin polarons. 
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